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ABSTRACT 

The interfacial stability of rotating core-annular flows is investigated. The linear and 
nonlinear effects are considered for the case when the annular region is very thin. Both 
asymptotic and numerical methods are used to solve the flow in the core and film regions 
which are coupled by a difference in viscosity and density. The long-time behaviour of 
the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a 
modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to 
three dimensions. The flows considered are applicable to a wide array of physical problems 
where liquid films are used to lubricate higher or lower viscosity core fluids, for which a 
concentric arrangement is desired. Linearized solutions show that the effects of density 
and viscosity stratification are crucial to the stability of the interface. Rotation generally 
destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces 
tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow 
finite amplitude helically travelling waves to exist when the fluids have different viscosities. 
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for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, 
VA 23681-0001. This work was also supported by SERC. 
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1 Introduction 


In 1967 Yih [20] studied the linearized long-wave stability problem for the parallel flow 
of two superposed fluids trapped between infinite flat plates. He showed analytically 
that Poiseuille and Plane Couette flows of viscously stratified fluids may be unstable to 
interfacial disturbances. Energy is transferred from the moving boundary or the applied 
pressure gradient to the wavy motion of the interface. Yih also found that the volume 
ratio of the two liquids is a crucial factor in determining the flow stability. Moreover, 
when the less viscous liquid layer is thin, the flow is found to be stable to disturbances 
with streamwise wavelengths much longer than the plate separation. 

Since this novel work there has been much interest in viscously stratified flows and 
their wide-spread practical applications. An example of this is the study of water waves 
generated by wind as considered by Blennerhassett [2]. He used Yih’s long- wave ap- 
proach to determine the growth of small amplitude surface waves for fully developed 
laminar flow of air over water, and then extended this to include the weakly nonlinear 
effects and determine the wave amplitude equation. The numerical results suggest that 
the nonlinear effects reduce the growth rate and lead to equilibrium amplitude waves. 

Another practical application of this theory arises from the difficulties encountered 
in transporting highly viscous liquids, such as crude oils. These difficulties can be 
overcome by heating, or diluting the fluid, such methods, however are often considered 
expensive and impractical, and a desirable alternative is to lubricate the flow using 
another immiscible fluid. Joseph, Renardy and Renardy [12] investigated the validity of 
the ‘viscous dissipation principle’ and concluded that a less viscous lubricating film has 
a tendency to encapsulate a heavy liquid core. 

Flows of concentric, immiscible fluids inside pipes are known as core-annular flows, 
hereinafter referred to as CAF. The most simple arrangement is when one liquid occupies 
the central core region and is surrounded by another filling the annulus. Our interest 
lies with such an arrangement inside a cylindrical pipe, known as a two phase CAF. 
When stable, the two fluids form concentric regions aligned parallel to the pipe wall. In 
practice such an ideal arrangement is not realised purely by chance and a question of 
some importance is to determine the circumstances under which this geometry is altered. 
For bicomponent flows of oil and water, (or some similar aqueous solution), in cylindrical 
CAFs practical observations include the formation of slugs of oil and water, wavy core 
flows and emulsions of both fluids. 

The linear stability of fluid-fluid interfaces has been studied comprehensively by a 
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number of authors. Hickox [7] considered cylindrical CAF in which the less viscous fluid 
occupied the core region and concluded that the arrangement was always unstable to 
long waves. Joseph, Renardy and Renardy [12] studied the inverse problem, that is, a 
CAF with a more viscous core, they too confirmed Yih’s findings (for a planar geometry), 
namely that the flow is stable to axisymmetric long wave disturbances. The effects of 
surface tension and gravity were neglected in this work, but for planar geometries the 
destabilizing influence of density stratification can be overcome by the stabilizing effects 
of surface tension and viscosity differences, provided that the depth of the less viscous 
fluid is made small enough, see Renardy [19]. Hooper and Boyd [8] indicated that it is 
short waves which are stabilized the most. 

The effect of surface tension on CAF is more complicated. Preziosi, Chen and Joseph 
[18] undertook a detailed linear stability analysis for an oil core flow lubricated by water 
in a two phase CAF, their findings are largely in agreement with the experimental work 
of Charles, Govier and Hodgson [3]. They conclude that surface tension stabilizes short 
waves and destabilizes long waves. The generalised capillary instability of long waves 
below the lower branch of the neutral curve is thought to lead to oil slugs and bubbles 
forming in the lubricating water layer. Preziosi et al. also investigated the influence of 
density stratification in the absence of gravity, they show that it has opposite effects at 
the upper and lower branches of the neutral curve. Increasing the density of the fluid in 
the annulus relative to that of the core destabilizes disturbances near the upper branch 
and to a lesser extent, stabilizes disturbances near the lower branch. 

Hu and Joseph [9] studied three arrangements of CAF: (i) An oil core with water 
annulus; (ii) A water core and oil annulus; and (iii) a three phase CAF with a water 
layer trapped between oil regions. They show that for the second case the flow is always 
unstable, but only weakly so when the oil film is thin. They also note that in all cases, 
when the Reynolds number is small it is the capillary instability which is most dominant. 

In much of the work mentioned above the influence of gravity is neglected. Both 
Chen, Bai and Joseph [4] and Georgiou, Maldarelli, Papageorgiou and Rumschitzki 
[6] considered its effect when the pipe is aligned vertically, so as not to violate the 
axisymmetry of the problem. The former authors found that gravity had a destabilizing 
effect at all Reynolds numbers with either a light or heavy lubricant. In particular 
they conclude that for slow flow, a heavy lubricant should be used when the pressure 
gradient acts in the same sense as the gravitational force, a lighter one when the flow 
acts upwards. Georgiou et al. [6] made a full analytical analysis of the linear stability of 


2 


CAF in which the depth of the outer layer is asymptotically small. They make several 
important conclusions, namely that the inertia of the film does not affect the stability 
at leading order and that for moderate surface tension it is the viscosity stratification 
which is dominant. There work was found to be in excellent agreement with the above 
authors when the depth of the film becomes larger. 

Miesen, Beijnon, Duijvestijn, Oliemans and Verheggen [14] found excellent quali- 
tative agreement between the theoretical and experimental study of interfacial waves 
in CAF. They explain that the absence of good quantitative results, especially when 
the wave amplitude was observed to be large, is most likely due to the ommission of 
nonlinear effects. 

In all this work so far we have discussed the case of a fixed pipe, and it has been 
found that the axisymmetric modes are the most dangerous. The effects of rotation 
on the stability of CAF are of practical importance, an example of this is in the oil 
industry where fluids are pumped through boreholes containing rotating drills. The 
study of a single rotating fluid has of course been undertaken by many authors, it is 
widely accepted that the most unstable disturbances are the axisymmetric modes. Hu 
and Joseph [10] performed a linear stability analysis for a two phase rotating CAF. They 
conclude that rotation stabilized the axisymmetric (n = 0) modes and destabilized the 
non-axisymmetric (n = 1) modes which were found to be the most dangerous. When the 
more dense fluid is in the annulus, there exists a critical rotational speed, above which 
the centripetal acceleration stabilizes the flow. With the lighter fluid in the annulus the 
flow is always unstable. 

In this work we study the weakly nonlinear interfacial stability of a rotating two 
phase CAF, for both axisymmetric and non-axisymmetric disturbances to the interface 
between the core and annular regions. The procedure used here closely follows that of 
Papageorgiou et al. [15] who studied the fixed pipe axisymmetric problem. In particu- 
lar we consider the case when the annular region is very thin compared with the pipe 
radius, so that an asymptotic expansion can be made in terms of a single small param- 
eter. We aim to analyse the effect of rotation on the perturbation to the concentric 
basic flow when the viscosities of the core and film layers are different. We also allow 
for a small difference in the densities of the two fluids, but the effects of gravity are 
neglected. Frenkel, Babchin, Livich, Shlang and Sivashinsky [5] considered the simpli- 
fied axisymmetric case when the viscosities of the film and core were equal so that the 
dynamics of the core did not affect the nonlinear stability problem. They obtained the 
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Kuramoto-Sivashinsky (KS) equation which governs the evolution of the interface, it can 
be written in normalized form as: 

Vt + VVz + Vzz + Mzzzz - 0 , 

where t denotes time and 0 the axial distance along the pipe. 

The KS equation arises in many physical situations such as flame propagation and thin 
film flow down an inclined plane. Extensive numerical studies on the route to chaos taken 
by solutions of the KS equation on a spatially periodic domain have been undertaken 
by Papageorgiou and Smyrlis [16] and [17], They observe many interesting characteris- 
tics shown by this equation, such as, steady modal attractors, periodic attractors and 
chaotic oscillations. The accuracy of their numerical experiments have identified win- 
dows, (ranges in the parameter A), where the time periodic solutions undergo a complete 
sequence of up to thirteen period doubling bifurcations. This has enabled the authors 
to compute Feigenbaum’s universal constant to a high degree of accuracy. 

Papageorgiou, Maldarelli and Rumschitzki [15] included the dynamics of the core 
which are coupled to the film flow by the viscosity difference. They found that a necessary 
constraint for the inclusion of capillary effects was that eJ ~ R e where J is a non- 
dimensional surface tension parameter, R e the Reynolds number based on the core speed, 
and e the ratio of the film depth to the radius of the core. They satisfied this by 
considering two physical regimes: 

A: slow core flow with moderate surface tension, 

B: moderate core flow with large surface tension. 

The resulting evolution equation is a modified form of the KS equation, it includes 
an additional term which represents the effect of the core dynamics. For case (A) this 
term has a purely dispersive effect because the slow flow gives rise to a Stokes flow 
approximation in which the inertial effects of the core are neglected. When either the 
more or less viscous fluid occupies the film, they find that this dissipative term tends 
to organise the otherwise chaotic motion of the interface into travelling waves. For the 
case of much larger surface tension (B), they find that viscous stratification also has a 
linearly stabilizing influence when the core is the more viscous of the two fluids, for a 
less viscous core it is linearly destabilizing. Never the less the imaginary contribution 
still has this same ability to produce travelling waves. 

When rotation is introduced the method used by Papageorgiou et al. can be adapted 
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to find the solution to the non-axisymmetric problem. We solve the core problem numer- 
ically as well as analytically for two specific cases, namely, a slow Stokes flow without 
rotation and secondly when rotation is asymptotically large. Having taken a thin film 
limit, a lubrication approximation is valid in the film layer and the solution in that 
region follows. The method by which the core dynamics are coupled with the motion of 
the film is through the equations of tangential stress, namely the radial derivatives of 
the axial and azimuthal flow. 

The capillary effects are retained if the speed of the core satisfies a similar constraint 
to the one discussed above. As we shall see later however, this is altered somewhat by 
rotation, the cenripetal pinch effect on the interface requires an even slower core flow 
to realise capillary effects. In other words, for a given Reynolds number, (that is, for a 
given axial speed), increasing rotation reduces the effect of surface tension. 

A solution is found for the disturbed flow in both the core and film regions and 
the evolution equation is derived. For the axisymmetric case, associated with rotating 
cylindrical CAF it is not surprisingly found that the nonlinear evolution equation which 
includes capillary and viscous effects is also a modified form of the Kuramoto-Sivashinsky 
equation. The rotation alters the additional linear viscous coupling term found by 
Papageorgiou et al. and its effect is discussed later. The non-axisymmetric problem 
extends the evolution equation to a novel three dimensional form of the modified KS 
equation. 

As rotation increases an interfacial boundary layer is formed inside the core layer. 
The analysis used here requires the thickness of the film to remain smaller than the 
depth of this layer, as a result we have an upper bound for the rotation of the cylinder. 

The layout of this work is as follows: in §2 we derive the non-dimensional basic flow 
for a rotating cylinder containing two concentric fluids; in §3 we take the limit of small 
film thickness and make a lubrication type approximation in the film. The core problem 
is then solved both numerically and analytically. The solution in the film is found and 
the evolution equation derived. In §4 we present the methods used to solve the linear 
and nonlinear problems and discuss the numerical results. Finally, in section §5 we draw 
some conclusions. 
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2 Geometry and Basic Flow 


Two incompressible, viscous fluids are contained in an infinitely long circular cylinder. 
The fluids have constant, but different densities and kinematic viscosities, and they 
are immiscible. We use cylindrical polar coordinates (r*, 0 *, 2 *) to denote the radial, 
azimuthal and axial directions respectively, so that the cylinder has radius R 2 and its 
axis lies along r* = 0. The interface between the two fluids is at a constant radial 
distance R\, so that they occupy two distinct, concentric regions: 

CORE REGION : Fluid of density p c and kinematic viscosity p c occupies 0 < r* < Ri, 
FILM REGION : Fluid of density pj and kinematic viscosity pj occupies R x < r* < R 2 . 


The cylinder rotates with angular velocity u and there is a prescribed, constant axial 
pressure gradient — F, (such that F > 0) but the effect of gravity is neglected. We seek an 
exact solution of the Navier Stokes equations of the form U c j(r*) = (0, V , W ) c j and 
pressure ~P* c j(r m ), where the subscripts c, / denote a core and film quantity respectively. 
There is a discontinuity in the film and core pressure fields due to the surface tension cr 
between the two fluids, (see Lamb [13]), hence: 

T'{r' = ft) -W = R, ) = |-. 

The basic flow is: 


v c (n 

= uir* , 

Vj(n 

= a;r*, 

W c (r-) 

= y~(R\ 



= F ( R l 

Ap s \ 1 


F 

4 Pj 


We non-dimensionalise the above quantities using a velocity scale Wo, the primary 
flow evaluated at the axis of the cylinder, 


W 0 = W* c (r* = 0) 


FR\ 


+ 


(r\ - ri) . 


4 Pc 4 pj 

The interface radius R x is used as a radial and axial length scale, whilst ^ and p c W£ 
are the appropriate scales for the time and pressure respectively. We also define 


R F — 


H'o R,\ p c 

Pc 


ft = 


u)R\R e 

Wo 
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N n _PJ n R2 

m = — , p - — , a 

Pc 


Pc Pc 

Here R e is the Reynolds number and fi is a rotation parameter. Using non-dimensional 
polar coordinates defined by (i r*,0*,z *) = (R x r,0, Riz), the primary flow may be written 
in the form 


— Hr 
V«(r) = -jp 

— Hr 
V/(r) = 


mr 


Wc{r) = 1 - T7 7, 

w a 2 + m - 1 

W/fr) 


a 2 — r 2 


(2.1a) 

(2.1b) 

(2.1c) 

(2. 1 d) 


a 2 + m — 1 

Consider now, a more general fluid velocity ([/, V, VU) C ,/ and pressure P c j which are also 
functions of the azimuthal coordinate 9. The non-dimensional governing equations in 
the core are: 

VU e V 2 _ 1 U 2V e 


u l + uu r + ^ + wu,-^ = -R + ?-?)• 
y 1 + t /K + 12i + w, + ^ = -T + i^-^ + Tr)- 


tt ? £ T UW r + 


and in the film: 


r 


+ W* 


- -r, + ^(v’iy), 


Rf/fl V 2 

Ut + UUr + — + WU Z 

r r 


U 2V* 


v 2 t/ - — - 


Wi 

y t + UV T + — + WV Z + 


_Pr_ 

p p R 
UV P e . m 


j »2 ^2 


= — - + 


vy + t/vy + 


r 


+ HTy 


pr pR e 

P z m 


+ 


,2 p 


P + P R e 


(v 2 w),, 


(2.2a) 

(2.2b) 

(2.2c) 

(2.3a) 

(2.3b) 

(2.3c) 


where the Laplacian operator V 2 (.) = ((.)rr + + (■ )zz + ) an d the subscripts 

denote partial differentiation. 

The equation of continuity must be satisfied by both the core and film fluids and takes 
the form 


(Ur) r + V$ + rW g = 0. 


(2.4) 


The film flow must satisfy no slip at the cylinder wall, (r — a), and in the core the 
velocity and pressure must remain finite at the axis. We now allow the interface to be 
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expressed more generally by r = S(0,z,t). Here we must satisfy continuity of velocity 
so that, 

Uc = Uf, V c = Vf, W c = W f . (2.5) 

For the continuity of tangential and normal stress we first resolve into directions parallel 
and perpendicular to the interface, we use the jump notation [.] = (,) c — (,)y, evaluated 
at r = S to obtain: 


+ (1 + Sl)F, - 


[R Z K + F Z \ = 0, 
F z S z S e I 


inhere in the core region: 

* ■ ¥(-' 

* = ¥(#' 
F, = JhL l(u.A 


• _ Sefl _ <7 F 1 _ 

1 T g 0 Z 1 2 


V Ue 
s + s 


= 0, 
J 


(V.n) , 


2 . A V Ue We 

PFe - -g (U + V 0 )j + Vr ~ - + -j - S a -± - S Z V Z 


(u, + W r - J- + k) + S z R e P — 2S z W^j , 


and the stress components in the film region are obtained from the above by replacing 
R e by R e m~ l . In the normal stress equation there is a discontinuity due to surface 
tension as discussed previously, here we have defined a non-dimensional surface tension 
parameter J = g ^’/ c which has been used by many authors including Hu and Joseph [9], 
Preziosi et al. [18] and Papageorgiou et ah [15]. The unit vector normal to the interface 
surface is denoted n and its divergenge is, 

„ a _ ( 1 r S„ , 2SI , SI S„si S„Sl , 2 S.S,S,.\ , i_3 
V.n - + — +___ 55- + — y3—J|n| , 

l n l = (l + + S. 2 ) ■ 


In addition to the above, we have the kinematic condition at the interface: 


-(r-S(0,z,t)) = 0 


u = St + WS z -{ 50, at r = 5(0, z, t). 


The full nonlinear equations (2.2a) to (2.3c) and boundary conditions (2.5) to (2.8) gov- 
ern the incompressible core-anular flow exactly, but in general they require a numerical 


8 



solution. Hu and Joseph [9] used a finite element method to solve the linearized energy 
problem for a general core to film ratio. In this work however, we restrict our attention 
to the case when the film layer is very thin compared with the radius of the core. This 
arises when the core fluid is separated from the pipe wall by a thin lubricating fluid 
layer. The method we use here essentially follows that of Papageorgiou et al. [15] by 
finding asymptotic solutions to the core and film fluids and hence we derive the non- 
linear amplitude equation which describes the evolution of the interface. In this work 
we investigate the effect of rotation and also allow for non-axisymmetric disturbances to 
the basic flow. 

3 Derivation of the Evolution Equation 

3.1 Thin Film Limit 

We take the radius of the cylinder to be given byr = a = l + e where e = + 0, so 

that the width of the film is O(t) and a lubrication type approximation is valid in that 
region. Suppose the general interface r — S is a small deformation of the previously 
constant position r = 1, hence we write r = S(0,z,t) = 1 + 6i](0,z,t), where ^ < f. 
Although r] is an unknown function, the interface conditions can be expanded in a Taylor 
series about the unperturbed position r = 1 . 

In the film we define a new variable y, such that: 

r = 1 + 7 — — ^nr, where y — 0(1) as e — > 0. (3.9) 

( 1 -p efl 3 J 

In the core we define: 

r = — — — where x = 0(1) as t — > 0. (3.10) 

1 T 0 3 

The motivation for the scaling factors involving will be made clear later when we 
consider the limit of large rotation. 

The basic velocity in the core and film is then perturbed by small amounts denoted 
(t/, V)W) c j and the surface tension across the interface induces a jump in pressure 
between the core and film fluids, denoted P c and Pj. We substitute these quantities 
into the governing equations of motion and interface conditions and discard those terms 
which are asymptotically small as we take the thin film limit e — ► 0. As we shall see, 
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weakly nonlinear effects do not enter the leading order momentum equations in the 
film and core but in fact do so through the interface conditions leading to a nonlinear 
evolution equation which governs the behaviour of the interface with time. 

Since we are interested in the effects of viscosity and density stratification between 
the two fluids, the core and film dynamics must be coupled together. This is achieved 
by the continuity of tangential stress, in particular, we balance the radial derivatives 
of both the axial and azimuthal film velocities with their corresponding counterparts in 
the core region. We define $ = (1 + Oa) and T = (1 + cOa), so that ~ ~ and 

dr ~ ~^c'Sy * n core anc ^ respectively. If we assume that the radial gradients of 
the azimuthal and axial velocity components are comparable then we have: 


dW c dWj_ 

dr dr 


W c 


4>e 


W f , 


dK dVj_ 

dr dr 



Evaluating the primary flow, (see equations (2.1a-d)) at the interface r = ^ = \ + Srj, 
we see that the azimuthal velocity is continuous but that there is a discontinuity of size 
O(S) in the axial flow, hence we require, W c ~ 5 . The core pressure scaling is obtained 
from the Navier Stokes Equations (2.2a-c), and the radial velocity size is found from 
the equation of continuity, (2.4). The film quantities all follow similarly. The perturbed 
basic flow is expanded in the form: 


n 6 

e 

rb 

kO 

U c — + . . . , 
$ 

u = u f n + 

Or 

Or 6<i>e 

— + V C 8 + ..., 

v =rP v ‘ r + 


— — 6$e 

W c + WJ + . . . , 

Wj = Wj + W f — + 

„ 0 V „ s$ 2 

„ p0 2 r 2 „ 6 $ T 

Fj+ 2Rl +P ‘R c +■■■ 

’ P - 2RI + P ‘ cR, + 


(3.11a) 
(3.11b) 
(3.11c) 
(3.1 Id) 


Where for convenience we have now dropped the tilde notation and Pj — J R~ 2 is the 
unperturbed pressure jump across the interface. In the view of these expressions, the 
kinematic condition yields, 


<f>e 2 




0 $e6Vj 


— Uj = TJ t + WfTJ z + —Wj T) z + I — + 


This can, however be simplified by redefining the time as t = 
formations 

z —> z - W j 0 t, 


o 


/Ip 


x = $ (1 + 6rj ) 
j and making the trans- 


it) 



where W / 0 is the basic axial velocity evaluated at the unperturbed interface accurate 
to 0(e ) , see [15]. This is equivalent to defining a coordinate system which is moving 
helically down the core, that is, one which rotates with speed 0R~ l and moves with 
speed W jq down the axis of the cylinder. We then obtain the asymptotic size of the 
interface displacement 6 in terms of the film width, so that 8 = and the kinematic 
condition reduces to: 


Uf{y =1) = T) T - —rjr) z . (3.12) 

m 

As mentioned above, it is the surface tension which produces a jump in pressure across 
the interface, equation (2.8) indicates that the film pressure is of the size In the 
light of equation (3. lid), we therefore require: 


Je 

Re 


~ 4>T ~ (l + 0 3 + efi 3 ) . 


(3.13) 


By isolating this capillary stability poblem, we are effectively restricting our attention 
to the phyisical regimes for which the speed of the core, characterised by the Reynolds 
number, and the surface tension J, satisfy the above criterion. When the cylinder is 
fixed, or rotating with negligible speed, 1, we see that (3.13) reduces to = 0(1), 
as found by Papageorgiou et al. [15]. To satisfy this they considered two distinct cases: 


A: a slowly moving core with moderate surface tension: R e ~ e and J ~ 1, 
B: a quicker core flow with very large surface tension: R e ~ 1 and J ~ e -1 . 


3.2 The Solution in the Core 


After perturbing the basic flow in the core, by the quantities given by (3.11a-d), the 
resulting leading order equations which govern the core flow are a linear system of 
partial differential equations which depend on radial distance x , azimuthal angle 0 and 
the axial coordinate z. The equations are written in terms of the new azimuthal, axial 
and time variable, thus for example we find that the inertial terms for the radial core 
momentum, transform into 


U t + 


_ 

Ug + WU Z 

X 


SUr + (W c 


Wjo) u z + 




u z + ... 



Ug, 


In order to solve the equations which govern the flow in the core we firstly employ a 
double Fourier transform to convert the partial differential system into one dependent 
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only on x and k, n the axial and azimuthal wavenumbers. 

We define: 

^ f + 7T y+oo 

U, V,w t p)= / {(U, V, W, P ) exp [-i ( kz + 710)]} dz dO. 

The governing equations in the core then reduce to 


<r v 2 f/- — - 


U 2 


_ $ 4 P' = 1 1 - 


X 2 X 2 


x‘ 

<j>2 


- 2 $ny, 


<h 3 vv 


y 2inU\ m$ 4 P 


+ ^J-— = ^-pj-W + 21W. 


$ 4 (v 2 w) - 2xR e U + ik&P = ^ ik$ 2 R e W, 

u m$y . , 

U H 1 1- tfcW = 0, 

X X 

where (.)' = ^ and V 2 (.) = ( ~ d ^ ^ ~ £(•))• 


(3.14a) 

(3.14b) 

(3.14c) 

(3.14d) 


These must be solved subject to the interface conditions: 

U = 0 = y, W=—(m-l), at x = 4>. 


W = — (m- 1), 

m 


(3.15) 


Here rj has ben defined by 

/ + 7 T r + oo 

/ {77(0, z, r)exp [— * (fcz + n 0 )]} dz d 0 . 

-7T — OO 

At the axis x = 0 we require the velocity and pressure to remain finite. 

As we shall see later, the focus of our interest in the core solution lies with the 
tangential stress components evaluated at the interface x = $: 


— (m — 1) N(0,z) = 


m 


* 2 (tv = )„ + (c«)« + ^ 2 (v) +( c/ ‘)»’l ( 3 - 16 > 

\ X / xQ 


In Fourier Space this can be written: 

C+ir f+oo 


„ ^ r-r7r y+oo 

N(n,k)= / / {N(0,z)exp[—i(kz + n8)]} dz dd. 

J — 1 r J — oo 

In general the solution of equations (3.14a-d) is a numerical task, and is discussed in 
section (3.3). We can, however, find analytical solutions for certain specific cases. 
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3.2.1 Slow Core Flow Without Rotation 


Consider the stability problem with moderate surface tension J ~ 1 as e — *■ 0, with the 
cylinder either fixed, or rotating with negligible speed < 1, hence $ ~ 1. Condition 
(3.13) indicates that we are interested in a slow core flow, namely R e ~ t. This corre- 
sponds to the equivalent axisymmetric problem (A) solved by Papageorgiou et al. [15]. 
The slow flow means that the inertial terms are neglected, so that in vector form the 
core equations reduce to V 2 (U c ) = VP, taking the divergence yields Stokes’ Equation 
V 2 P = 0. 


If we now make the transformations in Fourier space, Uj = U + iV and Vj — U — iV 
and define s = kx, then equations (3.14a) to (3.14d) reduce to 

nP 


P' - — - * v »+. (Pr) = 0 , 


p' + — - tv;., (V T ) = 0, 
iP-kVl(w) = 0 , 
VI (P) = 0 , 

U' T 4 Vi 4 I’-t . -l . ??: 4 H ZppL 4 2,4V = 0, 


where (.)' = ^ (•) = ® 


ds 


S 


(•))■ 


The interface conditions become 

U T = 0 = V T , W = — (m — 1) , at « = Jfc, 

m 

Hence the solution to the core flow follows as: 

U T (k,n,x) = Afz/ n (*)- 

V T (k,n,x) = A |z/ n (s) - ZzL^zliflj , 

W T (k,n,x) = iA |x/ n+1 (s) - | + frMf), 

m 2iT]kI n -i I n +\ 

[2(2 + n)/ n _i/ n+ i/ 7l — kl^l n - 1 — kI^I n +\ + 2kI%+ l I n ^i] 

7 n (s) is an n th order modified Bessel function of s and I n is its value at the point 5 = k. 
The function of stress components is then found to be: 

ifrv n 2in I n [kl n 2 7iI n +\I n kI n I n +\\ 

} ” - kP n In- 1 - klU^i + 2(2 + n)I n In+lln-l] 
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+ 

+ 


lk ^n+1 [kl n I n -1 2(?7 2^I n —\I n j r \ kI n I n +\ 

[2 */’ +1 / n -, - kP n I n ^ - */’/ n+1 + 2(2 + n)J n J n+1 J^] 

ikn. 


(3.17) 


For the axisymmetric case, n = 0, we have an identical result to the kernel 2 iNs found 
by Papageorgiou et al. [15]: 

M(0 k) 

3.2.2 Core Problem With Asymptotically Large Rotation 

Let us return briefly to the perturbed Navier Stokes Equations in the core, the reason 
why we rescaled the core radial variable in terms of U will now become clear. When 
rotation is large, 1, the Coriolis terms in the core momentum equations become 

dominant, and an analytical solution follows. To satisfy the inhomogeneous matching 
conditions at the interface, a non-trivial solution must be obtained near r = 1, we do 
this by increasing the size of the viscous terms in this region. 

We require the following asymptotic balances: 


Q,V ~ R e Pr , Vrr ~ ReP& ~ QU, 


W rr ~ R e P z , W ~ U r ~ V. 


Clearly then, ^ ~ S7 3 } this is the motivation behind the choice of x, as given by equation 
(3.10). We now let <f> = x — 1 — so that the interface of the two fluids corresponds to 
<j> — 0. The effect of this rescaling is to split the core fluid layer in two; the inner region 
contains stationary fluid and is matched asymptotically as <j> — > —00 to the outer core 
region which has thickness 0(fi~5). It is clear, however, that the following analysis is 
only applicable to very thin films in which e<C 1, since we require e <C Equation 
(3.13) indicates that for capillary stablility analysis we require R e ~ JefU 3 . Using 
the same method of Fourier transformation as in the previous section we find that the 
resulting core momentum equations are: 


2V = P\ 

V" = 2U + inP ; 
W" = ikP , 

U f = ~inV - ikW \ 


V {vi) = 4 k 2 V. 


which must be solved subject to the interface conditions: 


W = o) = o, Pv = 0) = , ^ ( ! ( ,t = 0) , v-w - 0) - 4,-t?(m - 1) 


2 k 2 


m 
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The solution is matched to the trivial one in the inner core layer as <f> -> -oo. 


Jl+iy/3) \ ( , (l — i\/3) 

V = A exp | (2k)* - J -</>\+Bexpl (2 k)*± L <t> 


+ C exp { (2 Jb)4) k > °> 

, (i -iVz) 


— A* exp —(2 k) 1 




2 <^] + 5*exp |-(2fc)s 2 




+ C*exp (2k) 3 (fj k < 0, 


^ = 
B = 
C = 


— — [3A; 2 — n 2 -f i(n 2 + & 2 )\/3 - 4mfc| , 

_l t2\ L v ' J 


-3A: 2 + n 2 + i(n 2 + fc 2 )"s/3 + 4mA;] , 


4m\/3(^ 2 + k 2 ) 

2 rj(m — 1) 
4m\/3(n 2 + k 2 ) 
—if) (m — 1) 

2m ’ 


where * denotes complex conjugate. 


The stress component function for k > 0 is then found to be: 

VY n (2A;) 3 [n + 3ik] 

N(n, k) = , 


(3.18) 


for non-positve k we have iV(n, —k) = N*(n , -\-k). 


3.3 Numerical Solution 

The solution of the general core problem (3.14a-d) is a numerical task. The conditions 
at the interface are given by equation (3.15), but since the equations are linear, and 
all other conditions are homogeneous, we normalise and solve the related system which 
satisfies: 


U(x = $) = 0 = V(x = $), W(x = $) = l. (3.19) 

Before proceeding further we must consider more closely the boundary conditions at the 
axis of the core x = 0. These have received much attention by many authors including, 
Batchelor and Gill [1] and Preziosi et al. [18], and so we briefly summarise their findings 
here. The physical quantities (U,V,W, P) c should remain finite at the origin and must 
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also be single valued there, that is, independent of 0 (see Joseph [11]). 

n = 0 : U( 0) = V(0) = W'(0) = 0, 

n = 1 : U( 0) + 1/(0) = W(0) = P(0) = 0, (3.20) 

n > 2 : U( 0) = 1>(0) = W(0) = P(0) - 0. 

The sixth order system of equations exhibit a singularity at the origin and admit three 

irregular solutions there which are to be discarded. We find the three remaining regular, 
independent power series solutions near x = 0 which satisfy the above conditions for the 
relevant azimuthal wavenumber n. Recurrence relations for these series are included for 
completeness in Appendix A. 

In the usual way, pressure is eliminated from equations (3.14a-d) so that they may 
be solved by means of a fourth order Runge Kutta scheme. The first N terms in the 
series are evaluated at x = £ -C 1 so that a numerical solution can be found by marching 
in x with step length h. A linear combination of the three solutions is then found to 
satisfy the interface conditions (3.19) at i = $. The size of the starting point £, the 
number of terms in the series N and the step length h are chosen so that this is achieved 
to within a required tolerance. The function of tangential stress components N(n, k ) is 
thus found to high degree of accuracy. The results are discussed in the next chapter. 


3.4 The Solution in the Film 


The film layer has thickness 0(e) and lubrication theory is appropriate there. We note 
that the film is thin enough to satisfy fi <C e -3 , as discussed previously, and that the 
equations are not directly affected by rotation at leading order since the Coriolis terms 
are negligible, and T ~ 1 . The solution which satisfies no slip at the cylinder wall 
(y = 0) is as follows: 


Pf = 


(3.21a) 

Wj = 

2 

l-P li+ yA(e,z,r), 

(3.21b) 

Vf - 

£-P li+ yB(e,Z, T), 

(3.21c) 

u s = 

+ m + j(At + B e >. 

(3.21d) 


The functions A and B which are independent of y are found by substituting the core 
solutions into the equations of tangential stress continuity: 




-1 

m 



mi 

«J)2 ’ 


(Wf) y = — 
y m 


TO. + 


(Uc) t 


$2 
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evaluated at y = 1 and x = $. 


3.5 Derivation of the Evolution Equation 

Now that we have solved the velocity and pressure fields for both the core and film layers, 
we can derive the equation which governs the time evolution of the small disturbance 
to the fluid-fluid interface. The kinematic condition is given by equation (3.12) and the 
continuity of normal sress is 


P]{y — 1) — (v + V zz + Vee) + /?o(l — p)v • (3.22) 


Here we have defined the surface tension parameter to be Tq = an 0(1) constant as 

e — » 0. The term involving Ro — due the density stratification of the two fluids 

and it arises as a result of the Taylor expansion of the basic pressure about the constant 
interface position. Clearly, when the rotation f) is very large this term will dominate 
equation (3.22), this means that the density stratification has a crucial effect on the 
evolution of the interface. We therefore restrict ourselves to flows in which the density 
difference between the two fluids is small, so that the term (1 — p)T 0 will always be at 
most an order one quantity. As discussed previously, we have chosen to neglect the effect 
of gravity for this problem, and so the role of density stratification is not as significant 
as the difference between the viscosities of the two fluids. The effect of gravity has been 
solved for the linear problem by Chen et al. [4] and Georgiou, Maldarelli, Papageorgiou 
and Rumschitzki [6]. 

By equating expressions for the radial film velocity Uj given by equations (3.12) and 
(3.21d), evaluating them at the unperturbed interface radius y=l and substituting for 
A and B, we obtain the evolution equation in the form: 



2 _L T ° ( 

m, + ^ 0 ». 


m 




+ Vee + zzzz + yeeee + 2 ? -] zz gg ) 

(m — 1) rjN(9, z ) 


+ 


m 2 $ 2 


= 0 . 


(3.23) 


Making the fluids the same density and considering the axisymmertic case d$ = 0, with- 
out rotation, reduces this identically to the evolution equation found by Papageorgiou 
et al., namely the Kuramoto-Sivashinsky equation together with the additional term 
N(9, z) which is due to the dynamics of the core fluid, coupled by the viscosity stratifi- 
cation. 
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4 Solution of the Evolution Equation 


The time evolution of the interface between the core and film fluids is governed by 
equation (3.23), which may be written in the form: 


Vr VVz 

m 

r+oo /*+o o 


+ iL ^ v2 ('') + ^ vs (''+ v2 ’') + 


1 ) 

47r 2 m 2< I> 2 


r°° f + °° N (n, k) r 

J — oo j — CO J— 7T J —OO 


where: V 2 ?? = ( d zz + dee) rj. 


The second derivative terms are linearly destabilizing and compete with the stabilizing 
fourth derivates. It is this competition which keeps the nonlinear solutions from breaking 
up, although there is as yet no formal proof that the solutions must remain bounded 
for all time when we have dependence on this additional angular coordinate 9. Before 
considering the diverse array of solutions offered by this nonlinear partial differential 
equation, it is firstly instructive to examine the simpler problem where the nonlinear 
term is absent. We look for solutions which are periodic in both the axial and azimuthal 
directions. Clearly r] must be single valued in 6 so that r/(9,z,t) — r/ (6 + 2n, z, f), 
however the length of the period in z is not determined and is indeed an integral part 
of the dynamics of the solution of the full nonlinear equation. 

We normalise the axial coordinate so that solutions, 77, are also 27T periodic in z, we then 
look for modes of the form ij ~ exp ct + iX^kz + inO^ , the linearized dispersion relation 
is: 

c r = a ( Xk 2 + n 2 ) + (AA; 2 + n 2 ) (l — A k 2 — n 2 ) — 7 N T ^k,nj . 

The parameter a = (1 — p) RqTq' represents the effect of the density stratification and 
rotation of the fluids, whilst the surface tension term has been normalized to unity by 
rescaling the time derivative. The core and film dynamics are coupled by the difference 
in viscosity between the two fluid layers and this enters the dispersion relation through 
the term 7 N r = where the subscript r denotes the real part. The sign of 7 is 

made negative when the core fluid is more viscous than the film. 

In accord with the usual theory of temporal stability, c is complex valued and a 
positive real part indicates exponential linear growth and hence instability. When both 
the viscosities and the densities are continuous across the interface, (a = 7 = 0), we can 
see from Figure 1 that the axisymmetric, n = 0 mode has a positive growth rate and is 
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the most unstable disturbance, in fact the azimuthal modes become progressively more 
stable with increasing wavenumber n. When the cylinder rotates, the density ratio of 
the two fluids can alter the stability of the CAF. If the core fluid is more dense than the 
film, ( a > 0), the non-axisymmetric perturbations become destabilized, with the n = 1 
mode being most dangerous, see figure 2a. For the inverse arrangement, that is with the 
lighter fluid in the core then all disturbances are stabilized due to the centripetal effects 
as shown by figure 2b. This is in agreement with the conclusions of Hu and Joseph [9]. 

More importantly we wish to quantify the influence of the core fluid dynamics on the 
interfacial disturbance. For a given surface tension and rotation we use the numerical 
results obtained from section (3.3) to calculate the tangential stress at the interface and 
hence the term N{k , n). In addition to the numerical solution of the core problem we can 
use the analytical results of sections (3.2.1) and (3.2.2). In the former case we have a slow 
core flow and negligible rotation leading to a solution of the Stokes’ problem for which 
we obtained a term N(k,n) which is imaginary, and hence has a purely dispersive effect 
on the evolution. For the simplified axisymmetric case this reduces to the kernel 2iNs(k) 
found by Papageorgiou et al. [15]. Although the linearized solutions are unaltered by 
this term, it does play an important role in the nonlinear evolution. 

In the presence of larger rotation we obtain the much simpler kernel : 


jy ( t, n) = mll|±M 
N (t,n) = ( 2 *> i[n - 3 ’*], 


k > 0, 
k < 0. 


Figure 3 shows that when the core fluid is more viscous than the film, (m < 1 and 
hence 7 < 0), the interface is stabilized according to linear theory. When the core is 
less viscous, for example the model representing the flow of lubricated crudes through 
rotating drills, we see from Figure 4 that the first non-axisymmetric mode is unstable 
and the linear analysis predicts exponential growth of the interfacial disturbances. 

This clearly is an unphysical situation and we must turn to the nonlinear theory 
for an explanation of such effects. Practical observations of such bicomponent rotating 
flows indicate that stable traveling waves may exist as well as the chaotic mixing of the 
transported fluid and its less viscous lubricant. 
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4.1 Numerical Method 

For the solution of the full nonlinear evolution equation we follow the method used in 
the extensive numerical work done by Papageorgiou and Smyrlis [16] and [17] in their 
comprehensive investigations of the route to chaos taken by spatially periodic solutions 
of the Kuramoto-Sivashinsky equation. 

Ut 4* UU X + U X x + ^Uxxxx = o, 

l/(x,0) = U 0 (x), U{x + 2ir,t) = U(x,t). 

The numerical method used in this work essentially follows theirs and consequently we 
need only summarise it briefly and discuss its extension to include the additional terms 
which represent the azimuthal coordinate. As described above we look for solutions of 
the evolution equation which have period 2w in both the azimuthal and axial directions 
having made the transformation of the form z — > \*z. 

We used a pseudo-spectral method in which the spatial dependence of rj is represented 
in wave number space by employing double Fourier transforms. The initial condition 
used by Papageorgiou and Smyrlis was given by U 0 = — sin(x), here we considered two 
separate initial value problems in which rj is a prescribed function of z and 0 of the form: 

CO OO 

7 / = jr r) k ,n sin (kz) cos (nO) , (4.24) 

k=\ n= 0 

o° °° * 

rj = 12 {nkl™s(kz + n6) + sin (fcz + n0)} . (4.25) 

1 n=-oo 

A traditional Galerkin approximation was then used to produce the following truncated 
system of ODE’s for 1 < k < K and -N <n < N. These are solved in a time marching 
sequence using a Runge-Kutta scheme. 

= hk tn T]k <n T Fk tn (flit, n) 

The term linear Ajt, n arises from the Fourier transform of the second and fourth derivative 
terms together with the kernel, as discussed above. The nonlinear contribution rjr] z is 
found by manipulating the series expressions for 77 , these may then be expressed as a 
combination of the initial spatial profile for each of the two cases considered. For 77 of 
the from (4.24), the term F k<n can be written as: 

£ + £ + E 

*l +*2 = * ki+k 2 =k k!-k 2 =k 

71 1 *h n 2 — n 11 

1 71 1 —tl 2 \—n ni+n2=n 
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kiVkun i*Ik 2 ,n 2 sin (kz) cos ( nO ) . 


- E + E 

k 2 -ki=k \k } +k 2 \-k 

7T.i+7i 2 = 7i |ri i —n 2 |=n 


E 

k 2 — k\ = k 
|m— n 2 |=n > 


For (4.25) we have: 


Fk,n = cos (kz + nO) - jFl*l sin (kz + nO) 


pW 

^ kyTl 


7 »( 5 ) 

fc,n 


2 *,n v ’ / 4 

= V n (c) n (5) 4- V r„ (c) T 7 (,) -n (c) r, (s) 1 

2 ^ V/ci ,nj Vk 2 ,n 2 ' Z-/ [Vk 2 ,n 2 'lki ,n\ Vki ,n\ T fk 2 ,n 2 J 

*f & 2 =& ki~k 2 =k 

71 ] -j-Tl 2 = 71 71 j — 71 2 = 71 

= V fn (c) n (c) - n (s) 77 (s) 1+T2 fn (c) r, (c) 4- n {s) n {s) 

/ y [ ? /fc 1 ,ni 7fc 2l n 2 ,ni 7fc 2 f n 2 ] ' ^ [% 1 ,m'/fc 1 ,ni ' 7fci ,n i 7fc 2 


712 I 


fcl + fc 2 =fc 

71 j +71 2 =71 


*1-A: 2 =fc 

71 1 —71 2 =71 


Throughout the numerical solution we compute the energy, (or L 2 -norm), the Fourier 
modes and a phase-plane map, which is the energy of the system against its time deriva- 
tive. An adaptive time step was used so that the error was minimised and the solution 
remained stable. For a full discussion of these points see [16]. 


4.2 Numerical Results 

Let us consider first the solutions obtained when the densities and viscosities are the 
same for both the core and film fluids, so that p — 1, and m = 1. We used a prescribed 
function of the form (4.24), with initial conditions which contained both axisymmetric 
and non-axisymmetric modes. When A > 1 constant states are obtained, as A decreases 
below unity however an array of nonlinear phenomena are observed. The conclusion of 
the lengthy computational studies undertaken in this work is that the azimuthal modes 
decay rapidly and only the axisymmetric, n = 0, modes are significant. Consequently 
the results obtained here are closely related to those found by Papageorgiou and Smyrlis 
[16] and [17] for the solution of the standard Kuramoto-Sivashinsky equation. 


Parameter Range 

Description 

1 < A < oo 

Constant States 

0.248 < A < 1 

Fully Modal Steady Attractor 

0.0755 < A < 0.247 

Bimodal Steady Attractor 

0.034 < A < 0.075 

Fully Modal Steady Attractor 

A < 0.033 

Mostly Chaotic Oscillations 
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Table 1 . Summary of the Solutions with matched densities and viscosities. 

The results illustrated in table L are not to be considered as a full description of the 
solutions of this equation, to do this we would require many hundreds more numerical 
experiments. Moreover, we seek to indicate that the rough estimates to the window 
bounds shown above are consistent with those for the standard KS equation and that 
no new behaviour has been discovered. When A < 0.033 the solutions become chaotic, 
although a closer numerical analysis indicates that time-periodic behaviour may also 
exist for sub-windows within this region of parameter space. As A decreases, more 
accurate, and hence longer, numerical experiments must be made to define precisely the 
behaviour as the solution undergoes a transition towards full chaos. 

Similar numerical studies were undertaken using a solution of the form given by 
equation (4.25). In the absence of a viscousity differnce between the fluid layers we see 
similar solutions to those outlined above, in fact if the cosine behaviour is dropped the 
solutions are identical. 

Figure 5(a) shows a plot of energy against time for a particular case in the chaotic 
window, namely A = 0.02. We note that although the solution remains bounded for all 
time, the structure is indeed that of fully chaotic oscillations as the phase-plane diagram, 
as figure 5(b) indicates. 

By preserving the modes r)^ n and 77 ^ at successive time intervals, we may visualise 
the evolution of the interfacial disturbance. The chaotic solution as a function of periodic 
axial variable z, and fixed azimuthal coordinate 6 , is shown by figure 6. The solution is 
shown at successive time intervals, and these are plotted with a displacement of 1 unit. 

We now consider the evolution of the interface when the two fluid layers are strat- 
ifed by a difference in their kinematic viscosities, so that we introduce the kernel term 
proportional to N(d,z). This additional term given by the coupled dynamics of the 
core region can alter the stability of the interface considerably. As we have discussed 
previously, this term is found either numerically for a general core speed, surface tension 
and rotation or may be given by the analytical work in (3.2.1) and (3.2.2). 

In this work we present the results for a quickly rotating core-annular fluid with a 
relatively low coefficient of surface tension between the two fluids layers. In this case the 
kernel N (n, k ) is complex valued and its imaginary part is proportional to k The real 
part allows for linear growth or decay of non-axisymmetric disturbences, as discussed 
above. 
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We have already seen that in the absence of viscosity stratification the interface may 
become chaotic in time, an important question then is what effect does the this viscosity 
difference make when the full nonlinear evolution equation is solved? For a less viscous 
core layer, how does the nonlinear behaviour cope with the exponential linear growth of 
the non-axisymmetric disturbances to the interface? 

The results of our numerical investigations are that the viscousity difference organises 
the otherwise chaotic motion of the interface into periodic waves which travel in both 
the axial and azimuthal directions. After only a short amount of time the energy, (the 
integral of J] 2 over its domain), becomes constant. The solution is illustrated by figures 
7 and 8 which show the axial and azimuthal travelling waves respectively. In Figure 7 
we have displaced the graphs vertically by 2 units and we have included three periods 
so that the motion of the wave can be visualised more clearly. When analysing these 
solutions one must recall that this is in fact a wave superposed on a coordinate system 
which already moves helically down the core. Figure 8 clearly indicates that this new 
interfacial profile is indeed non-axisymmetric. Again we include three periods and have 
a vertical displacement of 15 units. 

Upon increasing the ratio of the core viscosity to that of the film fluid, the amplitude 
of these waves is increased, but does not grow exponentially in time as the results of the 
linear theory given here and that of Hu et al. [10] suggests. If the core is less viscous 
than the film fluid the linearly stabilizing influence could ultimately yield the trivial 
solution. 


5 Conclusions 

In this work we have studied the linear and nonlinear behaviour of a rotating two-phase 
core-annular flow and identified the important mechanisms which govern the instability 
of the fluid-fluid interface. 

When a rotating fluid encapsulates a more dense core, the centripetal forces lead 
to the growth of non-axisymmetric disturbances to the concentric arrangement. More 
importantly perhaps, a viscosity difference between the fluids also leads to an unsta- 
ble non-axisymmetric regime. The linear theory predicts that these three dimensional 
disturbances will grow exponentially in time. 

The numerical solutions of the weakly nonlinear evolution equation derived here 
indicate that in fact travelling waves are obtained. These waves are distortions of the 
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interface and move helically down the core. Their amplitude depends upon the difference 
in the viscosities of the two fluids. 
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A Power Series Coefficients 

In section 3.3 we discussed the numerical solution of the core problem and seek a power 
series solution to equations (3.14a-d) near the axis x = 0. We look for a solution of the 
form: 

oo CO 

v = £ u,r\ v = Y. vy, 

3 - 0 3=0 

oo oo 

^=x;wy, p = ev. 

j = 0 j = o 

In the light of the conditions given by (3.20) must consider each case separately. 


A.l Axisymmetric Modes n=0 

W 0 , Pq and V\ are arbitrary; U 0 , W x , Vo, P\, U 2 , V 2 , W 3 and P 3 are all zero. 

U u = -\ikW 0 , 
u 3 , = - l -ikw 2 , 


1 / 3 , = ^( 2 MJi + ikaR'Vi - k 2 V x ) , 

1 


IV 2 , — — (k 2 W 0 + ike? Pq + ikReWoj , 


P 2 , = 


2c x 4 


(8 c?U 3 - ikU x ~ k 2 U\ + 2afIVi) . 


For j = 4, 5, 6,. . . , 


Uj* = 
Vh = 


ikWj.i 

"(T+ry’ 

1 

a 2 (i 2 - 1) 


a 2 k 2 Vj - 2 + 2aQ.Uj- 2 + ika 2 R e V^ 2 - ikR e Vj- A 
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w h = -^\a 2 k 2 W j - 2 + ika 4 P j - 2 -2R e U j -3 

a 1 ) 1 

+ ikR e a 2 Wj-2 ~ ikR e Wj-4 , 

Pj-, = — 2 <* 3 f)Vj_i — ika 2 R e Uj-\ -\-ikR e Uj - 3 

Ot j 

+ i(i + 2)a 4 (/ i+1 -aW i _ 1 ]. 

A. 2 Non-Axisymmetric Modes 

For the first non-axisymmetric mode n — 1: Wj, Fi and Uq are arbitrary, Wo, Po, U u 
Vi, W 2 and F 2 are all zero. 

U 2 , = (k 2 Uo - 2iaZCtUo + ikR e a 2 U 0 + c* 6 ^ - 2 ikaW x ) , 

V 2 , = 3iU 2 -kW u 

W 3 , = ~ (a 2 k 2 Wi + ika 6 Pi - 2 R e Uo + ika 2 R e Wi) , 

U 3 , = 0, 

v 3 , = 0 . 

For subsequent modes n > 2: W n , P n and U n +i are arbitrary, f/o, Vo, Wo, F 0 ,. • ■ ,U n , 14, 
F n _i, W n -i are all zero. 

T/ _ i{n + 2)U n+ i kW n 

V n+ 1 , = 

n n 

W n+1 , = 0, 

Pn+l, = 0. 

For j = n + 2, n + 3, . . ., 

Uj, = r 2 2 ] [ a2 k 2 Uj - 2 - 2ifcWj_i + (j - 1) a 6 Fj_i 

[0 + 1 ) - ™ ] 

- 2a 3 f2Vj_2 + ikR e (< a 2 U } _ 2 - f/,_ 4 )] , 

= a , JJT r „ 2) + iiaVj., + .kR, (a 2 Wj-2 - W^) 

p j’ = ^j[in 2 + ((j 2 + 2j~n 2 )u^2))(ikR,( a %.,-U j . 3 ) 
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- 2a 3 QVj-i - 2 ika 4 Wj 

+ a 2 k 2 Uj -\ ) + ({j + 2) 2 — n 2 ) (nkaR e (V)_ 3 — a 2 Vj- i) + SlUj-i 
+ 2ia 2 ink 2 a 3 Vj-\ + ( j 2 + 2 j - n 2 ) ika 4 W , )] , 

/(j) = ((;' 4- 2) 2 - n 2 ) - (j 2 + 2j) (j 2 + 2> - n 2 ) + 2n 2 j + n 2 . 
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Fig. 3 : Linearised Evolution : Large Rotation Sc Less Viscous Core 
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Figure 8 Solution 77 in 
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